Stability and electronic structure of the (lxl) SrTiO3(110) polar surfaces by first 

principles calculations 
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The electronic and atomic structure of several (1 x 1) terminations of the (110) polar orientation 
of SrTiOa surface are systematically studied by first-principles calculations. The electronic structure 
of the two stoichiometric SrTiO- and 02-terminations are characterized by marked differences with 
respect to the bulk, as a consequence of the polarity compensation. In the former, the Fermi level 
is located at the bottom of the conduction band, while in the latter the formation of a peroxo bond 
between the two surface oxygens results in a small-gap insulating surface with states in the gap of the 
bulk projected band structure. We also consider three non stoichiometric terminations with TiO, 
Sr and O compositions, respectively, in the outermost atomic layer, which automatically allows 
the surface to be free from any macroscopic polarization. They are all insulating. The cleavage 
and surface energies of the five terminations are computed and compared, taking into account the 
influence of the chemical environment as a function of the relative richness in O and Sr. From our 
calculations it appears that some (110) faces can even compete with the TiC>2 and SrO terminations 
of the (100) cleavage surface: in particular, the (HO)-TiO termination is stable in Sr-poor conditions, 
the (llO)-Sr one in simultaneously O- and Sr-rich environments. The available experimental data 
are compared to the outcomes of our calculations and discussed. 

PACS numbers: 68.47.Gh, 31.15.Ar, 73.20.At, 68.35.Bs 
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I. INTRODUCTION 

Strontium titanate is a material that has many po- 
tential applications. Its high dielectric constant makes 
it a candidate to replace silicon dioxide in some nano- 
electronic devices. SrTiOa has been used in photo- 
electrolysis and as a substrate for the growth of high 
T c superconductors as well as of other thin oxide films 1 . 
Like other perovskite compounds, SrTiC>3 shows a rich 
physical behavior. It undergoes an anti-ferro-distortive 
(AFD) transition at T < 105K: the cubic structure turns 
into a tetragonal phase where the neighboring TiC>6 oc- 
tahedra are tilted with opposite angles around a [100] 
direction, doubling the unit cell. Moreover, SrTiC>3 is at 
the boundary of a ferroelectric (FE) transition 2 *^ but 
remains paraelectric for all temperatures^. 

While the (100) surfaces of SrTiC>3 have been exten- 
sively studied, both theoreticall y 6 ! 7 ! 8 , and experimen- 
talh y 9 ! 10 ! 11 , the (110) terminations are much less known. 
Such a scarcity is likely due to the polar character of the 
(110) orientation. The sequence of atomic layers of O2 
and SrTiO stoichiometry implies a monotonic raise of the 
microscopic electric field, which has to be compensated 
either through a modification of the surface composition - 
which leads to non-stoichiometric terminations - or by an 
anomalous filling of the surface electronic states - which 
must imply crucial variations of the electronic structure 
of the surfaces that should be in principle detectable by 
experiments. 

Nevertheless, {110} terminations of strontium titanate 
has been observed quite often, but the effort towards the 
precise characterization of their atomic-scale morphology 



and the detailed study of the corresponding electronic 
structures started only in the last decade. For instance, 
it is worth noting the investigations by atomic force 
microscopyiS, scanning tunnelling microscopy (STM), 
Auger spectroscopy and low energy electron diffrac- 
tion(LEED) measurements 1 ^, as well as ultraviolet pho- 
toemission and x-ray photoemission spectroscopies (UPS 
and XPS, respectively), coupled to LEED 14 . In the 
latter works the authors made a considerable effort to 
characterize the modifications of the electronic struc- 
ture that are connected to different preparation condi- 
tions of the surface through STM and scanning tun- 
nelling spectroscopy (STS) 1 ^, angle resolved photoemis- 
sioni^ and conduction measurements! 7 .. Not all these 
studies provide a unifying picture of the surface, which 
indeed shows a great sensitivity to the thermodynamic 
conditions, mainly temperature and oxygen partial pres- 
sure. 

Motivated by that experimental work, some theoreti- 
cal studies of SrTiO3(110) recently appeared. The op- 
timized surface geometry was obtained through inter- 
atomic forces that were derived either by an empirical 
shell- mo deli&, or by more accurate treatments that are 
based on the Hartree-Fock approximation, still describing 
the hamiltonian matrix on semi-empirical ground s 19 ! 20 . 
The approximate nature of those approaches allowed 
some trends to be deduced, but it did not permit in 
most cases to draw quantitative conclusions that could 
be compared to experiments. On the other hand, first- 
principle approaches that are based on the Density Func- 
tional Theory (DFT) have been widely applied to the 
study of oxide surfaces in recent years. More specifically 
to perovskite surfaces, such an approach has been showed 
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to be reliable and complementary to experiments^. It is 
the aim of the present paper to provide an application of 
state-of-the-art first-principle methods to the (110) po- 
lar termination of SrTi03, which is the first one to our 
knowledge. 

In particular, we focus on two issues: (i) Which is 
the mechanism at work to cancel the macroscopic po- 
larization? and (ii) Is the thermodynamic stability of 
the (110) termination comparable to that of the (100) 
cleavage face? Regarding (i), while in the wide-gap rock- 
salt compounds MgO and NaCl the non-stoichiometric 
reconstructions such as the octopolar one are clearly fa- 
vored over polarity compensation by an anomalous filling 
of surface states, it has been recently suggested for the 
less ionic ZnO(OOOl) and ZnO(OOOl) surfaces 2 ^ that 
the latter mechanism may be at work. The coexistence 
of the rather ionic Sr-0 bonds with the more covalent Ti- 
O ones, as well as the fairly large dielectric constant of 
SrTiC>3 might situate strontium titanate at the border- 
line between the two previous cases. Regarding (ii), the 
quite numerous observations of {110} (1 x 1) terminations 
suggest that these surfaces may be obtained in appro- 
priate thermodynamic conditions. However, the ternary 
nature of strontium titanate, as well as the lack of pre- 
cise experimental indications on the surface composition, 
make the first-principle modelling a rather delicate task. 
Therefore, in this work, we restrict ourselves to several 
terminations, stoichiometric or not, of (1 x 1) reconstruc- 
tions, completing a previous preliminary reporljS 3 .: the 
two stoichiometric O2 and SrTiO faces, and the three 
non-stoichiometric O, TiO and Sr terminations. Their 
atomic and electronic structures are computed and dis- 
cussed in Section III. Moreover, their relative stabilities 
are compared as a function of the chemical environment 
(see Sections IV and V). 



bohr and R c (p) — 1.75 bohr for the O, with s, p and d 
being the different channels. The Sr 4s and 4p semi- 
core electrons are treated as valence electrons in the self- 
consistent procedure, (which we refer to in the following 
as the small-core approximation) whereas the Ti 3s and 
3p semi-core states are frozen (large-core approximation) . 
Since the latter states have an appreciable superposition 
with the 4s and 3d radial orbitals of Ti, in the exchange 
and correlation functional we take into account the non- 
linear core corrections* between the radial density cor- 
responding to the frozen 3s and 3p Ti orbitals and the 
valence electron density. 

The previously described computational scheme is 
tested on the cubic phase of bulk SrTiC>3 . The lattice pa- 
rameter ao, the bulk modulus Bq and the cohesive energy 
E co h are very close to the experimental data (see Table 
P). The large-core approximation for Ti is also checked 
against the small-core one, and the former is adopted in 
all slab calculations. The gap between the valence and 
conduction band computed through the difference of the 
Kohn-Sham eigen-energies (2.1 eV) is smaller than the 
experimental value (3.2 cV) 30 , which is the common case 
for the DFT2i. 

TABLE I: Experimental and computed lattice parameter ao, 
bulk modulus Bo and cohesive energy E co h for the cubic 
SrTiC>3 phase. The experimental values of ao and E co h are 
quoted in Refp— and Bo is taken from Refp— . Two distinct 
calculations adopting the small-core and large-core approxi- 
mations for the Ti pseudopotential (see text) are also shown. 



Calc. (large-core) Calc. (small-core) Experiment 
00 (A) 3.951(+1.2%) 3.875(-0.8%) 3.903 
Bo (Gpa) 187(+2%) 192(+5%) 183 

E coh (eV) 34.06(+7%) 37.00(+17%) 31.7 



II. COMPUTATIONAL METHOD 

The calculations are carried out within the density 
functional theory (DFT)24. The exchange and correla- 
tion energy is treated via the local density approxima- 
tion (LDA) using the Perdew-Wang parametrization- 25 . 
The Kohn-Sham orbitals are expanded on a plane wave 
(PW) basis set and an energy cutoff of 30 Hartree (Ha) 
is employed. For all calculations, we use the ABINIT 
computer code2&. 

In conjunction with the use of a PW basis set, we 
adopt soft norm-conserving pseudopotentials that are 
generated following the Troullier-Martins scheme^ in 
the Kleinman-Bylander formSM, in order to avoid to take 
into account the inner atomic electrons in the self- 
consistent cycle explicitly. Along the pseudopotential 
generation process, atom reference states 3c? 1 4s 5 4p - 25 
for Ti, 4s 2 4p 6 5p 1 for Sr and 2s 2 2p 4 for O are used. The 
pseudization radii are i? c (s)=2.45 bohr, R c (p)=2.55 bohr 
and R c {d)=2.25 bohr for Ti, R c (s)2.0 bohr, R c {p)=2.5 
bohr and R c (d)=2.2 bohr for Sr and R c (s) = R c (d) = lA 



The calculated formation energies of bulk Ti and Sr 
from the isolated atomic species are: E^ = 4.58eV and 
e£ t = l.73eV, which differ by -6% and +2% from the 
experimental values, respectively. On the other hand, 
the O2 formation energy (E[, 2 = 7.58eV) is about 45% 
greater than the experimental value, like in previous LDA 
calculations^. 

The surfaces are described in the framework of the slab 
model. The calculations are carried out by sampling the 
irreducible Brillouin zone by a (4,4,2) Monkhorst-Pack 
mesbj 3 ^. By using a (5,5,2) mesh, the computed total 
energies decrease by less than 10 meV for the insulating 
terminations and 20 meV for the metallic ones. In those 
cases, when an effective Fermi-surface smearing is em- 
ployed to get a better convergence 3 ^, a careful procedure 
is used to extrapolate convergent energies in the limit 
T — > 0. The stable surface configurations are obtained 
by minimizing the Hcllmann-Feymann forces in the Born- 
Oppcnheimer approximation, generally starting from the 
geometry of the ideal unrelaxed surfaces. 
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In order to avoid spurious interaction between peri- 
odic slabs via dipole-dipole interactions^^, the compo- 
sitions of the adopted slabs are symmetric upon inversion 
along the surface normal. Well converged slab total ener- 
gies can be routinely obtained for typical vacuum widths 
ranging between 10 and 12 A. On the other hand, the de- 
termination of the minimal number of atomic layers that 
is allowed needs some care. For instance, a nine-layers 
slab is sufficient for the two non-stoichiometric (110)- 
Sr and (HO)-TiO terminations, whereas 11 layers are at 
least needed for the two (110)-O 2 and (HO)-SrTiO stoi- 
chiometric terminations. In the case of the O-terminated 
surfaces, a very careful choice of the bulk reference en- 
ergy is also needed, as detailed in the following. To 
this purpose, we made calculations up to 15-layers thick 
slabs. As a general rule, the slab is considered to be thick 
enough when the difference between the total energies of 
two slabs differing by a SrTi03 unit and the bulk total- 
energy is smaller than 1 mHA This condition cancels 
out the problem of nonconvergent surface energies^2il£ 
and makes it possible to obtain a fair bulk reference en- 
ergy for each termination. 

Even if our calculations are carried out at K, we im- 
pose symmetry constraints to our slabs to simulate the 
cubic phase that is stable at room temperature. A special 
case is represented by the non-stoichiometric (110)-O ter- 
mination, since the mirror symmetry along the [110] di- 
rection is lost. As it will be detailed in the section devoted 
to such a termination, the bulk energy should be com- 
puted with reference to an anti-ferro-distortive (AFD) 
phase that is simulated in a tetragonal supercell analo- 
gous to that used for the slab. Other authors gave details 
about the change of the AFD distortion with respect to 
the bulk structure^ 1 . Our computed total-energy differ- 
ences between the cubic and the AFD phases is about 2 
mHa. Neglecting such a tiny quantity may appreciably 
bias the numerical extrapolation of the surface energy 
when a linear scheme is used. To summarize, the un- 
certainty on the surface energies E c i within our slab ap- 
proach does not exceed 10 -2 J/m 2 , which is usually much 
smaller than the computed surface energy differences for 
the various terminations. 

We also focus on the total electron distribution in the 
various surface configurations. An estimate of the elec- 
tron sharing between the O anions and the Sr and Ti 
cations is obtained by means of the Bader's topological 
analysis of the electron density, which corresponds to a 
partition of the total charge in atomic bassins. Such a 
procedure is independent of the basis set useda^, and has 
been recently used^ 3 . and improved within the ABINIT 
package. The precision of the atomic charge integration 
is about 5 x 10 -3 electrons. 



III. ELECTRONIC AND ATOMIC STRUCTURE 

In the (110) orientation of SrTi03, a stacking sequence 
of atomic layers ...02-SrTiO-02-SrTiO... is provided. If 



we consider that the ionic charge of O, Ti and Sr are 
Qo = —2, Qti = +4 and Qsr = +2 respectively, then the 
O2 plane bears a formal charge Qo 2 — ~4 and the SrTiO 
plane QsrTiO = +4 per 2D unit cell. According to the 
criterion for polarity compensation 4 ^^, the formal sur- 
face charges of the various SrTiO3(110) faces have to be 
equal to half the bulk value, i.e., Q sur f = ±2. Therefore, 
we consider two main classes of (1 x 1) terminations: on 
one side, the so called stoichiometric terminations, since 
their compositions reflects the bulk stacking: the (110)- 
SrTiO and -O2 ones, for which an anomalous filling of 
surface states is expected. On the other side, we study 
the (HO)-TiO, -Sr and -O terminations, for which the sto- 
ichiometry changes can in principle provide the polarity 
compensation. For each termination, we describe in the 
following the detailed atomic and electronic structure. 



A. Stoichiometric terminations 

1. The SrTiO termination. 

In the stoichiometric (HO)-SrTiO termination (see 
Figure^), 6 Sr-0 bonds and 2 Ti-0 bonds are cut. Since 
the surface plane bears a formal charge QsrTio = +4 
there must be an anomalous filling of surface states, in or- 
der to cancel out the macroscopic component of the slab 
dipolc 44 . Figure |2 reports the computed band structure. 
One can note that the Fermi energy is located above the 
bottom of the conduction band. Some conduction-like 
states are thus filled, which bring the two additional elec- 
trons that are needed according to the electron counting 
rule for polarity compensation. 



o ° 




FIG. 1: Side view of the SrTiO slab cut along a [001] plane. 
Sr, Ti and O atoms are white, grey and red, respectively. 

These states are localized in the outermost layers and 
are genuine surface states. The density of the filled con- 
duction band state at the M point of the Brillouin zone 
that is displayed in figure |3 is roughly delocalized over 
three surface layers, which shows that the self-consistent 
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Lines in the Brillouin zone 



FIG. 2: Computed band structure for the SrTiO termination 
of SrTiOs(llO). Only the valence band top and the conduc- 
tion band bottom are drawn. 



charge redistribution is more complex than the model 
that is based upon the values of the formal ionic charges. 
A more thorough insight can be obtained by means of the 
Bader's topological charge analysis, which is reported in 
Table ITT1 The layer charge q^Tio an d <io 2 of the first and 
third SrTiO, and the second O2 outermost planes, respec- 
tively, are indeed strongly modified by the filling of the 
additional surface state and thus show marked changes 
from the computed charge for the inner bulk-like layers 
(q Bulk = ±2.50), which is roughly achieved only in the 
fourth layer. The polarity compensation criterion is ful- 
filled, since the sum of the charge of the four first layers is 
equal to -1.30, which is almost the half of the bulk layer 
charge q Bulk . 




SrTiO SrTiO SrTiO SrTiO SrTiO 

O. O, O, O, 



FIG. 3: Localisation of the compensatory state at the M 
(0.5;0.5;0.0) point of the Brillouin zone 

Finally, we can see in tabled that the difference be- 
tween the surface layer charge qg r Tio an( i the bulk layer 
charge q Bulk is mainly localized on the Titanium, which 
is confirmed by the direct visualization of the filled state 
at the conduction band bottom (in the right part of fig- 
ure 0}. However, an effective decrease of the electronic 
kinetic energy may be obtained through the derealiza- 
tion of this state over neighboring sites. Such a screening 
phenomenon, reminiscent to that found in Na-covered 
Ti02(f 40) surfaces^, thus affects the Sr and O charges, 
too. The charge of the outermost Ti is thus reduced, 



TABLE II: qo, qTi and qs r are the Bader's topological charge 
of the Oxygen, Titanium and Strontium respectively, for the 
SrTiO termination. I, II, III and IV give the index of the 
layer. Charge are given in electrons. For comparisons, bulk 
values are given at the bottom. 



layer 


Atomic charges 


Layer charge 


SrTiO 

2 
SrTiO 

2 


qo=-l-36; q^ = 1.74; q Sr =1.40 
q '=-1.38 
q" 7 =-1.30; qTi 7 =2.03; q|"=1.55 
q o V =-1.30 


qsrTio = l-78 
qo 2 =-2.76 

qsrTio=2-28 
q£=-2.60 


Bulk 


qgu' fc =-1.26; q ««** =2 .l8; q£" ifc = 1.58 


qt ,uU = 2 5Q 



which can be probed by XPS experiments through a sur- 
face Ti signal that should differ from the bulk Ti 4+ formal 
oxidation state to a great extent. 



Ti 




Ti 




FIG. 4: Cutting of the total valence density (left) and com- 
pensatory state density (right) perpendicularly to the [100] 
direction, in a Ti02 plane 

As far as the atomic structure is concerned, we can see 
in table lllll that large relaxations happen in the four out- 
ermost layers. For instance, the large rumpling between 
Oxygen and Strontium that is equal to 0.73 A in the sur- 
face layer, is still non negligible in the third layer (0.32 
A). The inter-planar distance between the outermost Sr- 
TiO layer and the O2 layer beneath contracts by 11 %, 
while that between the latter and the third SrTiO planes 
grows by +19%, showing a typical damped oscillation 
that makes the surface unlike layers closer two-by-two, 
which is accounted for by the general theorj^Siil. How- 
ever, the surface Ti-0 bondlengths cipi-o ar e only weakly 
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TABLE III: Relaxation and rumpling on the (HO)-SrTiO 
termination. The mean positions of the SrTiO layers are 
computed by averaging the normal coordinates of the cor- 
responding atoms. The interplanar distances are given in A 
and their relative variations in brackets. 



layer 


Relaxations and rumplings 


SrTiO 


Ti(-0.03) O(+0.38) Sr(-0.35) 


1 


1.25 (-11%) 


2 


no rumpling 


I 


1.66 (+19%) 


SrTiO 


Ti(+0.06) O(+0.13) Sr(-0.19) 


I 


1.28 (-8%) 


2 


no rumpling 


I 


1.49 (+6%) 


SrTiO 


Ti(+0.01) O(+0.04) Sr(-0.05) 


1 


1.38 (-1%) 


2 


no rumpling 



modified: for Ti in the topmost layer, dTi-o is reduced 
of -3%, while it is expanded by 5% in the third layer. 
Such a behavior is due to the combined relaxation and 
rumpling that mainly rotate the Ti-0 bond. At odds, 
the more ionic and weaker Sr-0 bonds are much more 
affected: dg r i_o"i is reduced by 13%, while doi-Sx m is 
expanded by 24% 



2. The O2 termination. 

The complementary stoichiometric (110)-O2 termina- 
tion is also expected to undergo deep modifications of 
the electronic structure because of the polarity compen- 
sation, which requires two electrons to be removed from 
the surface. As displayed in the left part of Figure^] the 
two surface oxygens move close to each other forming a 
bond 1.48 A long. Such a bond- length is usually a sig- 
nature of a peroxo group which is formally denoted 
as 0\~ in the ionic limit. It would imply that the sur- 
face charge per unit cell is Q SU rf = — 2 instead of -4 like 
in the bulk. Therefore, the formation of a peroxo group 
formally permits to fulfill the polarity compensation cri- 
terion. 

Less formally, one can see that the formation of a per- 
oxo bond is made possible by emptying the anti-bonding 
a* molecular orbital (see the right panel of Figure |SJ). 
This also provides an effective mechanism to open a small 
gap at the surface (see figure [SJ • The Fermi level rests 
above the two anti-bonding states 7r* and does not cross 
the band structure . The bulk gap value of SrTiC>3 
(around 2 eV) is found between the valence band and 
the conduction band (i.e. excluding the O2 tt* states). 
The two in-gap n* states are split since the two [110] and 
[001] crystallographic directions are not equivalent, and 
are localized at the surface layer as displayed in figure [7| 

Table IIVI yields the results of a Bader's topological 
analysis for this termination. The topological charges of 




FIG. 5: Left part: side view of the O 2 -SrTiO 3 (110) slab, cut 
perpendicular to the [001] crystallographic direction. Right 
part: schematic molecular diagram of an oxygen dimer. In 
the case of an isolated group, the two n* orbitals are 

completely filled. 
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FIG. 6: Computed band structure for the 02-terminated slab. 
Only the valence band top and the conduction band bottom 
are drawn. The ■k* states discussed in the text give rise to four 
bands (two for each termination) which are not completely 
degenerated due to finite-slab effects. 

the surface O (O 1 ) are sensitively reduced, consistently 
with the formation of the peroxo bond. We obtain a sur- 
face charge almost equal to half the bulk one. In this 
case, the charge modification is essentially restricted to 
the outermost O atoms, as a consequence of the covalent 
and localized nature of the peroxo bond. 

As far as the atomic structure is concerned, important 
modifications are present on the surface. The Oxygens 
tilt at the surface and reduce their angle with respect to 
the [110] direction by strongly increasing the bond- length 
with the Titanium on second layer (do I -Ti 11 ) by 19 %■ 
As a consequence, there is a big inward relaxation (-20% 
- see also table |Vj of the peroxo bond with respect to the 
second layer. On the other hand, d i_ Sl .ii is remarkably 
reduced (-13 %). 

We point out that a few surface configurations were 
obtained for the (110)-O2 termination through the min- 



r I 1 

1 11 i 

o I =^ '—- ' -— 

02 02 02 02 

z axis in atomic units 

FIG. 7: Localisation of the highest occupied state it* for the 
2 termination of SrTiO 3 (110). 



TABLE IV: Same caption as Tabic ITU for O2 termination 



layer 


Atomic charges 


Layer charge 


o 2 

SrTiO 

2 
SrTiO 


q o =-0.66 
q" =-1.20; qTi=2.13; q|£=1.55 

qo"=-1.21 
q v =-1.24; q^=2.19; q^=1.58 


qo 2 =-l-32 
q^ TiC .=2.48 
q^=-2.42 

ClSrTiO— 2.53 


Bulk 


qg u '*=-1.26; q^"'*=2.18; q# r u, *=1.58 





imization procedure. The one reported in Fig|S]has the 
lowest total energy. Another configuration with the per- 
oxo group on top of the Titanium is found, which is very 
slightly higher in energy (5.10~ 3 J/m 2 ) but still within 
the intrinsic precision of our calculations. These two con- 
figurations have very similar band structures and atomic 
charges, and can be considered as almost degenerated. 
It is worth noting that when starting the geometry opti- 
mization from the ideal, unrelaxed (110)-O2 termination, 
a third local minimum without peroxo groups at the sur- 
face is found. With respect to the peroxo configurations, 
it has an open-shell structure and a higher surface energy 
of about 0.7 J/m 2 . 



B. Non-stoichiometric terminations 

The (HO)-TiO, -O and -Sr terminations may be ideally 
obtained from the stoichiometric surface by adsorption 
or desorption of Sr and O atoms. They bear a formal 
charge Qno — +2, Qo = ^2 and Qg r = +2 in the ionic 
limit, respectively, which corresponds to half the bulk 
layer charge Qbuik = ±4. Therefore, they are all com- 
pensated and no anomalous filling of surface states is in 
principle necessary. In order to demonstrate this assump- 
tion, the computed density of states of those three non- 
stoichiometric terminations are drawn in Fig|SJ They 
are all insulating, with electronic structures qualitatively 
similar to the bulk. 
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TABLE V: Same caption as Table ITTTI for O2 termination 



Layer 


Relaxations and rumplings 


O2 


no rumpling 


1 


1.12 (-20%) 


SrTiO 


Ti(-0.05) O(+0.12) Sr(-0.07) 




1 52 f+9%1 


O2 


no rumpling 


1 


1.32 (-6%) 


SrTiO 


Ti(-0.02) O(+0.03) Sr(-O.Ol) 


1 


1.45 (-4%) 


O2 


no rumpling 


1 


1.40 (0%) 


SrTiO 


no rumpling 



TiO 




1 p 


1 ^ 
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If- 


Sr 
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A A* 


\ 


it 



-30 -20 -To o 



energy in eV 

FIG. 8: Density of states of the TiO, Sr and O terminations. 
The occupied states are filled in gray and Ef is the Fermi 
level. 



1. The Sr and TiO terminations. 

The (llO)-Sr termination is obtained from the O2 one 
by adsorbing a Sr atom per unit cell whereas the (110)- 
TiO one is obtained from the (HO)-SrTiO termination 
by removing a row of surface Sr along the [100] direction. 
These two terminations are displayed in figure El 

These two (110) terminations are representative of 
more open surfaces. The coordination numbers of the 
surface atoms are reduced with respect to the stoichio- 
metric (110) (1 x 1) terminations as well as to the non- 
polar (100) surfaces. Let us imagine to cleave a SrTiC>3 
slab along different orientations: two complementary ter- 
minations are obtained, for which a certain number of 
bonds are missing. Such cut bonds are: 4 Sr-0 and 1 
Ti-0 for the non-polar (100) orientation; 2 Ti-0 and 6 
Sr-0 bonds for the (HO)-SrTiO and (110)-O 2 termina- 
tions, as well as for the (110)-O termination, which is self- 
complementary; 2 Ti-0 and 8 Sr-0 bonds for the (110)- 
TiO and (llO)-Sr terminations. Moreover, at variance 
with the Ti0 2 and SrO (100) terminations and SrTiO, 
2 and O (110) ones, the (llO)-Sr and the (HO)-TiO ter- 
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TABLE VII: Same caption as Table El for the (HO)-TiO 
termination. For the (llO)-Sr termination, the topological 
charges of the surface atoms are very similar to the bulk ones, 
from which they differ by 7% at most. 



layer 


Atomic charges 


Layer charge 


TiO 

Oa 
SrTiO 
2 


q^=-1.20; q^=2.04 

q£>"=-1.20; q^f =2.21; qlf =1.58 
q^=-1.25 


qTio=0.84 
qo,=-2-22 
q^i io =2.59 
qg:=-2.50 


Bulk 


qg u,fc =-1.26; q^ u,fc =2.18; qf r ui *=1.58 


q W *=-2.50 



FIG. 9: Left (right) panel: side view of the TiO (Sr) slab, cut 
along a [001] plane. 



minations are the only ones having an under coordinated 
atom even in the third layer starting from the surface (a 
Sr-0 bond is cut). Therefore, it is not surprising that the 
remaining atoms on the (llO)-Sr and (HO)-TiO undergo 
large relaxations (see table IVIf) . like the inward relax- 
ation of the surface Sr on the Sr termination (-43%) and 
the strong rumpling of the Ti and O atoms on the TiO 
termination (0.48A). Despite this under-coordination, 
relaxations are more quickly damped than for the stoi- 
chiometric terminations. For instance, the inter-planar 
relaxation between the third and the fourth layer is re- 
duced by a factor two with respect to the stoichiometric 
(HO)-SrTiO and (110)-O 2 terminations. 



TABLE VI: Same as TableED for the (HO)-TiO and (HO)-Sr 
terminations. 



Layer 


Relaxations (TiO) 


Relaxations (Sr) 


TiO/Sr 


Ti(-0.24) O(+0.24) 


no rumpling 


I 


1.51 (+8%) 


0.79 (-43%) 


O2 


no rumpling 


no rumpling 


I 


1.43 (+2%) 


1.53 (+10%) 


SrTiO 


Ti(+.08) O(-.07) Sr(-.Ol) 


Ti(+.06) O(-.02) Sr(-.04) 


I 


1.32 (-5%) 


1.35 (-4%) 


O2 


no rumpling 


no rumpling 


I 


1.40 (0%) 


1.41 (+1%) 


SrTiO 


no rumpling 


no rumpling 



Since the (llO)-Sr and the (HO)-TiO terminations are 
intrinsically compensated by stoichiometry, we may guess 
that the atomic surface charges differ little from the bulk, 
and that such an effect is mainly be due to the reduc- 
tion of the surface atom coordination numbers. In fact, 
the Bader's topological analysis (see Table WT^i confirms 
the polarity compensation of the two terminations: the 
charge modification affects essentially the three outer- 
most layers and their sum on these layers is equal to 
+1.21, which is almost half of the bulk value. Even 
if these changes are less intense than on the stoichio- 
metric terminations, they extend rather deeply into the 
slabs, which confirms the correlation with the presence of 



under-coordinated atoms and the interplay between the 
atomic relaxations and the electron redistribution. 



2. The O termination. 

Such a termination may be ideally obtained in two dif- 
ferent ways: by removing an Oxygen atom from a O2 ter- 
mination or by adding an Oxygen atom on a SrTiO ter- 
mination. The presence of a single O atom on the surface 
breaks the mirror symmetry along the [110] direction, at 
odds with all (1 x 1) terminations that were previously 
considered. Consequently, among all the (110) termi- 
nations, the convergence of the atomic structure and the 
surface energy of this one is the most delicate. We can see 
in figure HTH that the broken surface symmetry permits a 
distortion and its propagation into the slab. It consists in 
an alternating rotation of the octahedra along a [100] di- 
rection, which is known as an anti-ferro-distortive (AFD) 
instability. We point out that such a abnormal deep prop- 
agation of the AFD instability from the surface into the 
bulk is related to the fact that the low-symmetry phase is 
stable at T=0 K, the temperature for which the calcula- 
tions are performed. This propagation may be hindered 
by considering a 15 layer slab with the five inner layers 
frozen. Even in that case, the surface atomic structure is 
practically undistinguishable from that obtained on top 
of an AFD bulk. 




FIG. 10: Side view of the (110)-O slab. 



8 



TABLE VIII: Relaxation and rumpling on the (110)-O ter- 
mination. The mean positions of the SrTiO layers are com- 
puted by averaging the normal coordinates of the correspond- 
ing atoms. The inter-planar distances are given in A and their 
relative variations in brackets. <5Ti0 2 corresponds to the Ti-0 
lateral relaxation (the projected interatomic distance along 
the [110] direction). 





R plfi"vat i on ti anrl T'liTnnlino' 


^Ti02 









I 


0.94 (-33%) 




SrTiO 


Ti(-0.09) O(+0.06) Sr(+0.03) 


0.391 


1 


1.42 (+2%) 




o 2 


O(±0.33) 




1 


1.38 (-1%)- 




SrTiO 


Ti(+0.05) O(+0.02) Sr(-0.08) 


0.087 


1 


1.38 (-1%) 




02 


O(±0.23) 




1 


1.376 (-1%) 




SrTiO 


Ti(-0.03) 0(0.00) Sr(-0.03) 


0.057 


1 


1.40 (0%) 




o 2 


0(± 0.25) 





The surface O relaxes inwards (-33%) and laterally tra- 
in 



ward the two Strontium underneath (see table IVILTfl 
the inner part of the slab, no large inter-planar relaxation 
happens. However, a quite strong rumpling affects the in- 
ner O2 layers, which corresponds to the AFD distortion 
and is not damped as a function of the slab thickness. 
Moreover, there are large lateral atomic displacements 
within the SrTiO layers that are denoted with <5Ti0 2 ( see 
tablc lVliTll and decrease quite rapidly going into the bulk. 
As a consequence of this distortion, the O-Ti-0 aligne- 
ment is broken along the [001] direction. 
The topological charges of the surface oxygen's are very 
similar to the bulk ones, from which they differ by 10% 
at most. 



A. Cleavage energy 

When a stoichiometric SrTiC>3 slab is ideally cut 
and the two parts are put apart, two complementary 
surface terminations are created. In our case, they 
are the (HO)-SrTiO and the (110)-C>2 terminations on 
one hand, the (HO)-TiO and the (llO)-Sr ones on the 
other. The (110)-O termination can be considered as 
self-complementary. The respective cleavage energies 
-E , d° 2+SrTl0 ' > and E^J l0+Sl ^ can be obtained from the to- 
tal energies computed for the symmetric slabs through 
the following equations: 



E 



(TiO+Sr) 
cl 



— ( F ni 

2>S* 



TiO 
b 



E 



(0 2 +SrTiO) 
•I 



——(F U2 

\ slab 



^slab 



pSrTiO 
^slab 



nEbulk) 

- nEbulk) 



(1) 
(2) 



where E l slab the total-energy of the symmetric slab with 
the i termination, n the total number of bulk formula 
units in the two slabs, S the surface area and Ebuik the 
bulk energy per formula unit in the cubic structure. 
Since the (110)-O termination is self-complementary 

and shows an AFD distortion, its cleavage energy E^f + °' 
is calculated by using a unit factor instead of one half and 
a different reference energy E^f, which reads: 



E 



0+0 

cl 



S 



(E° ab - nE, 



AFD\ 
bulk ) 



(3) 



IV. THERMODYNAMIC STABILITY 



As previously pointed out in Section II, a different choice 
for the reference energy would result in a ill-defined cleav- 
age energy. This is exemplified in a recent calculation, in 
which an apparently diverging cleavage energy was ob- 
tained for the (110)-O terminationi& by using the cu- 
bic, undistorted reference energy instead of E^® . Our 
numerical results are summarized in table IIXI The O 
termination has the lowest cleavage energy (2.54 J/m 2 ) 
whereas the O2 and SrTiO the highest one (6.52 J/m 2 ). 
We also note that the three non-stoichiometric termina- 
tions have lower cleavage energies than the two stoichio- 
metric ones. Consequently, the polarity compensation 
that is achieved through the modification of the surface 
stoichiometry seems to be more effective than that by 
the anomalous filling of the surface states, as far as the 
energetics is concerned. 



Because of the different chemical nature of the various 
(1 x 1) (110) terminations of SrTiOa, their stability must 
be discussed as a function of the actual surface compo- 
sition 7 ii2ii2i^£. To this purpose, we make use of two dis- 
tinct physical quantities: the first one corresponds to the 
energy that is necessary to split the crystal in two parts 
and create two complementary terminations, which we 
refer to as the cleavage energy E c [. The second one is 
the surface grand-potential fi* , which is a measure of the 
excess energy of a semi-infinite crystal exposing a termi- 
nation with a given composition i in contact with matter 
reservoirs. 



TABLE IX: The cleavage energies E ct in J/m 2 . 

O + O TiO + Sr SrTiO +Th 
E c i 2.54 3.86 6.52 



B. The surface grand-potential. 

In order to distinguish the contribution of each termi- 
nation to the cleavage energy, we compute its surface 
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grand potential, which implies a contact with matter 
reservoirs. Many authors have recently used this method 
with success, as well for binary^S*^ as for ternary com- 
pounds**^. We introduce the chemical potential /XTi, fJ-sr 
and no of the Ti, Sr and O atomic species, respectively. 
The surface Grand Potential per area unit fl l of the i 
termination reads^: 

& = ^ IKlab ~ ^TiMTi - JVsrMSr - N »o] (4) 

with TVri, Nst and No the number of Ti, Sr and O 
atoms in the slab and the half factor corresponds to the 
surface Grand potential by termination. The chemical 
potential ^SrTi0 3 of a condensed and stoichiometric phase 
of strontium titanate is written as a sum of three terms 
representing the chemical potential of each species within 
the crystal: 

MSrTiOa = MSr + MTi + 3/-i (5) 

Since the surface is in equilibrium with the bulk SrTiOs, 
we have MsrTi0 3 == -Ebuifc- If we replace equation [S] in 
equation 0] we can eliminate the //Ti and /zsrTi0 3 vari- 
ables in the surface Grand Potential and obtain: 

= ^Kia6-^Ti^fe-M0(iV0-3iV T1 )-M Sr (iV S r-iV Ti )] 

(6) 

Relying upon equation^ one can deduce, for each ter- 
mination, the range of the accessible values of f2* if the 
minimum and maximum values of the O and Sr chemical 
potentials are known (see appendix). 

If we introduce the variation of the chemical potentials 
with respect to those computed for the reference phases 



E mol 

(Afi = Mo t~ and Afj, Sr = nsr-E^ lk , respectively) 

in equation we obtain: 



^ = & - —[Allo(N - 3N Ti ) ~ A[i Sl (N Sr - N Ti )} 
with fa = -^[E\ lab ~ N Ti E bu i k - 

-^-{No ~ 37V Ti ) - E*? k (N Sl - N Ti )] (7) 

fa measures the stability of the surface with respect to 
bulk SrTiOs, gaseous molecular oxygen and metallic Sr. 
In Table Ixl our ab initio results for fa, are listed. 

TABLE X: The surface energy fa as defined in Eq.|7]is given 
for each termination i in J/m 2 . 



i O TiO Sr SrTiO Q 2 
fa 1.27 6.88 -3.01 4.78 1.73 



For two complementary terminations (TiO-Sr and 
SrTiO-02) the sum of their surface Grand potential is 
independent of the chemical potential and corresponds 
to their cleavage energy (see Table llXll . 

The derivation of the upper and lower bounds of A/xo 
and A/isr is detailed in appendix A. Within the allowed 
region, we show in figure ITT1 (left panel) the (1 x 1) (110) 
terminations having the lowest surface grand potentials, 
which provides the stability diagram of SrTiO^llO) (1 x 
1) surfaces in an O and Sr external environment. 



An„(eV) k 




' O TiO(llO) -6 -5 -4 -3 -2 -1 

O poor 



Oxygen chemical potential A \l in eV 

FIG. 11: Stability diagram of the of the (lxl) SrTi03 (110) surface. The actual most stable termination is represented in the 
left panel as a function of the excess O and Sr chemical potentials A/i (vertical) and Afisr (horizontal). In the right panel, the 
surface Grand potentials are represented as a function of Afio (for the particular value of the Sr chemical potential Afisr ~ 
eV). 
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First of all, according to FiglTTl (left panel), our cal- 
culations predict that only four out of the five possible 
terminations may be obtained. Indeed, the (110)-O2 ter- 
mination cannot be stabilized, even in very O-rich chem- 
ical environments. The Sr termination is the most stable 
one in O- and Sr-rich environments, as its complementary 
TiO face is in O- and Sr-poor conditions. The O termi- 
nation shows a stability domain in moderate O- and Sr- 
environment. Finally, the stoichiometric and open-shell 
(HO)-SrTiO termination happens to be stable in a small 
domain corresponding to O-poor and Sr-rich conditions. 
In order to dicuss the existence of this small domain (at 
least within the theory) with respect to the precision of 
our calculations, we point out that (see the right panel 
of Fig. 11): 

• At the O-poor and simultaneously Sr-rich zone 
boundary (i.e.: A/io=-5.46 eV and A/igr= eV) 
the difference between Sr and SrTiO surface Grand 
Potentials is equal to 0.13 J/m 2 , which is higher 
than the estimated precision of our calculations, 
which is of the order of 0.01 J/m 2 . 

• Since the area of the (AfiQ, AfJ-s r ) domain is 
slightly underestimated by our first principle cal- 
culations (see appendix), the consideration of the 
experimental boundaries would enlarge the stabil- 
ity domain of the SrTiO termination. 

V. DISCUSSION 

Among the questions that we raise in the introduction, 
in the following we discuss three issues that are especially 
noteworthy. 

A. A stable polar stoichiometric termination: 
(HO)-SrTiO. 

We have shown in the previous section that the (lxl) 
stoichiometric (HO)-SrTiO termination can be stabilized 
in a O-poor and Sr-rich environment. Even if its stability 
domain comes out to be rather small, this is the second 
case (to our knowledge) together with the Zn and O ter- 
minations of ZnO(0001)2i of a stable stoichiometric po- 
lar oxide surface with an open-shell electronic structure. 
However, a recent STM investigation 51 questioned the 
proposed stoichiometric morphologies for ZnO(OOOl) and 
provides different, non stoichiometric, structural models. 

The a polarity compensation through anomalous fill- 
ing of surface states is expected on the basis of theoreti- 
cal arguments, but it is actually not often encountered^. 
Indeed, the cleavage energies of stoichiometric polar ter- 
minations in rock-salt structures, such as the (1 x 1) un- 
reconstructed (111) faces of MgO, are generally much 



higher than those of reconstructed, non-stoichiometric 
polar terminations. As a consequence, the surface grand 
potentials of these polar stoichiometric terminations re- 
sult higher than those of the non stoichiometric ones even 
when the chemical potential dependent terms represent 
a negative contribution. In this respect, the peculiar be- 
havior of SrTi03 may be due to the presence of Ti-0 
covalent bonds, in conjunction with a not too large fun- 
damental gap. Indeed, the energy increase that is due to 
the anomalous filling of surface states needed for polarity 
compensation may be effectively lessen by atomic relax- 
ation and electronic screening effects. The latter ones 
show up through a non negligible charge transfer affect- 
ing some surface layers as supported by the analysis of 
the topological Bader's charge of the SrTiO termination 
( Table Hill : all Ti and O charges belonging to the three 
topmost layers are modified. Such a screening mechanism 
is less costly than a drastic charge reduction on only one 
or two surface atoms, which is essentially the case for 
MgO(lll) (1 x 1)22, 



B. Comparison with the non polar SrTiOs(lOO) 
terminations. 

The cleavage (100) orientation of SrTiOa is non-polar 
and displays two different stoichiometric terminations: 
the (100)-TiO 2 and the (100)-SrO. They are usually ex- 
pected to be more stable than any polar face such as 
the (110) terminations. Indeed, within the same theo- 
retical and computational framework, we calculated the 
cleavage energy of the (1 x 1) (100)-TiO2 and -SrO sur- 
faces, which is not very much lower than those of the non- 
stoichiometric (110) terminations. Therefore, it is worth 
comparing the thermodynamic stability of the simulated 
(1 x 1) (110) terminations with the two (1 x 1) (100)-TiO 2 
and (100)-SrO faces. In figurelT^lthe domains of stability 
of the two (100) terminations and the five (110) termi- 
nations previously reported are gathered. Even if in the 
most common conditions, corresponding to moderate Sr 
and O chemical potentials, the (100) faces are favored, 
two (110) distinct (1 x 1) terminations are predicted to 
be stable - the (HO)-TiO in Sr-poor environments and 
the (llO)-Sr in O-rich and Sr-rich conditions - the first of 
which shows a quite wide domain of thermodynamic sta- 
bility. From the theoretical point of view, the (HO)-TiO 
and the (llO)-Sr terminations can be thought as exposing 
strongly relaxed {100} nano- facets (see Fig0J. It is there- 
fore not completely surprising that their surface energies 
can become comparable to those of flat {100} terraces, 
if the formation energy of corners and edges is not too 
high. This is consistent with the experimental evidence, 
which shows that {110} orientations can be quite easily 
obtained for SrTi03. 
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A|i (eV) 



-16.38 -9.23 -7.19 -2.58 




SrO(lOO) 

| TiO2(100) 



poor 

FIG. 12: Phase diagram: comparison between the different 
terminations of the (lxl) SrTiOa (110) SrTiO, 2 , TiO, Sr 
and O terminations and SrTiOa (100) Ti02 and SrO termi- 
nations in an Oxygen and Strontium external environment. 



C. Comparison with experimental results. 

Using the first-principles stability diagrams displayed 
above, we can now discuss and propose a tentative ex- 
planation for the experimental measurements. As an- 
ticipated in the introduction, the SrTiO3(110) surface 
was produced and characterized by several groups, show- 
ing a great sensitivity to the preparation conditions. In 
particular, if the SrTiO3(110) is slightly annealed in ul- 
trahigh vacuum (UHV), at 800°C i& or at 960°C dur- 
ing 2 hoursi£, the surface exhibits a (lxl) LEED pat- 
tern. In contrast, a (n x m) periodicity is observed 
when the SrTiO3(110) surface is annealed at tempera- 
tures T > 1000°C. Therefore, one may guess that (lxl) 
SrTiC>3 (110) surfaces represent local minima of the pos- 
sibly complicated free-energy landscape in a quite wide 
temperature range. Preliminary calculations on larger 
cell reconstructions^ show that some (lxl) phases may 
remain stable with respect to other reconstructed sur- 
faces. In the following, we mainly discuss the experiments 
performed on these (lxl) phases. 

As far as the composition of the surface is concerned, 
the authors of Ref. [l3| get evidence through Auger mea- 
surements of an increasing surface concentration of Sr 
at the high-temperature annealed surface. Therefore, 
they conclude for a lack of Strontium for the mildly an- 
nealed (lxl) surface, and propose a {100} micro-faceted 
(110)-TiO2 termination as the corresponding atomic- 
scale structural model. On the other hand, the STM 
study of Refi*3i shows on the (lxl) terminations the pres- 
ence of rather flat regions where tunnelling spectra have 
metallic character. These findings are consistent with 
typical Ti 3+ and Ti 2+ features in the XPS spectra, as 
well as a metallic surface state with a Ti 3d character as 



seen in UPS. On the basis of all those observations, and 
at odds with Ref. ITsl they propose a SrTiO termination 
for the unreconstructed SrTiO3(110) surface, which can 
account for its flatness and metallic character dominated 
by Ti-like occupied states. 

Relying on our calculations, we argue that these two 
apparently contradictory models are not necessarily in- 
compatible. On one hand, the UHV Auger measure- 
ments may have been carried out in Sr-poor environ- 
ments, whereas the UHV experimental conditions that 
have been used in Ref. 1 15l imply an O-poor environment. 
Keeping in mind the great sensitivity of the SrTiO3(110) 
surface to the actual thermodynamic conditions, it is not 
surprising that its atomic-scale structure may show big 
variations. According to our calculations (see Fig lllf) . 
the (HO)-TiO termination, which can be interpreted as 
a {100}-TiO2 micro- faceted surface in agreement with 
the model proposed by Brunen and Zegenhageni^, can 
be obtained in a wide domain corresponding to slightly 
Sr-poor and O-poor environments. If we assume that 
the temperature annealing is done in O-poor and, at the 
same time, rather Sr-rich conditions (whether one starts 
from the very beginning with a stoichiometric SrTiO ter- 
mination without desorbing Sr atoms or the latter ones 
migrate from the bulk to surface domains of TiO com- 
position remains an open question) a local (HO)-SrTiO 
termination may be obtained as proposed by Bando and 
coworkers^. Indeed, the metallic character of the (110)- 
SrTiO termination that is due to a surface state mainly 
of Ti character agrees with the measured UPS spectra. 
Moreover, as we have previously pointed out for this ter- 
mination, the anomalous filling of a Ti surface state that 
is needed for polarity compensation, should be associated 
to special Ti™ + features in XPS, with n < 4. 



VI. CONCLUSION 

The SrTiO 3 (110)(lxl) surface was studied in the 
framework of ab initio calculations for the first time. 
The number of surface terminations in this ternary 
compound is larger than in binary compounds, which 
permits a remarkable variety for mechanisms of polarity 
compensation. An anomalous filling of the surfaces 
states takes place at the SrTiO termination with an 
open-shell electronic structure and at the O2 termination 
with in-gap states. If non-stoichiometric O, TiO and Sr 
terminations are considered, we obtain for each of them 
an insulating bulk-like electronic structure. 

By calculating the surface grand potentials, we obtain 
four distinct (1 x 1) stable (110) terminations as a 
function of the Sr and O chemical potentials. A quite 
large domain is found for the three non-stoichiometric 
terminations, especially for the TiO one. This compe- 
tition between different terminations is also confirmed 
by the presence of a small domain of stability for the 
SrTiO termination in simultaneously Sr-rich and O-poor 
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